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^ . The mass function of clusters of galaxies is determined from 400 deg^ of early 

0\ . commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in 

^ . the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two indepen- 

O , dent selection methods: a Matched Filter and a red-sequence color magnitude 

Q . technique. The two methods yield consistent results. The cluster mass function 

^ ! is compared with large-scale cosmological simulations. We find a best-fit cluster 

^! normalization relation of ag^l^^'^ = 0.33 ± 0.03 (for 0.1 < 0.4), or equiv- 

2 ! alently as= (jf^)^'^. The amplitude of this relation is significantly lower than 

to '. the previous canonical value, implying that either Q^n is lower than previously 
expected {flm= 0.16 if (7s= 1) or ag is lower than expected (cr8= 0.7 if flm= 0.3) 
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as suggested by recent results. The shape of the cluster mass function partially 
breaks this classic degeneracy; we find best-fit parameters of ^m= 0.19 ±qq7 and 
as— 0.9 ±q;|. High values of 0,^ 0.4) and low ag (< 0.6) are excluded at > 

2(7. 

Subject headings: cosmology:observations-cosmology:theory-cosmological parameters- 
dark matter-galaxies:clusters:general- large-scale structure of universe 

1. Introduction 

The abundance of clusters of galaxies as a function of mass places one of the strongest 
constraints on the amplitude of mass fluctuations on 8 Mpc scale, cxg, and on the mass 
density parameter, flm. The present-day cluster mass function was the first observation to 
suggest that the standard ^m— 1 Cold Dark Matter (CDM) model has to be highly biased, 
with (T8~ 0.5 (i.e., a bias of ~ 2, since the galaxy fiuctuations amplitude is aggai ~ 1), in order 
to match the observed cluster abundance. The cluster mass function also showed that low- 
density CDM models fit the cluster data well with little or no bias (i.e., 1) (Bahcall & 
Cen 1992; White, Efstathiou, & Frenk 1993). The mass function constraint, frequently called 
"cluster normalization" because of its powerful constraint on the linear mass power spectrum 
amplitude (Xg, has provided the well known relation a^Qm^'^ = 0.5 ± 0.05; this result was 
obtained from observations of both the cluster mass function (Bahcall & Cen 1992, 1993) 
and from the cluster temperature function (Edge et al. 1990; Henry & Arnaud 1991; White, 
Efstathiou, & Prenk 1993; Kitayama & Suto 1996; Eke, Cole, & Prenk 1996; Viana & Liddle 
1996; Eke et al. 1998; Pen 1998; Markevitch 1998; Henry 2000). This relation implies that 
0.3 if (T8~ 0.9 - 1; this latter value of cig is suggested from other observations including 
cluster abundance evolution (Bahcall & Fan 1998; Donahue & Voit 2000, and references 
therein), the fiattening of the mass-to-light ratio on large scales (Bahcall, Lubin, & Dorman 
1995; Bahcall et al. 2000), and the SDSS and 2dF large scale structure observations (Szalay 
et al. 2002; Verde et al. 2002). Similar as-^m normalization relations have been recently 
obtained from weak lensing observations on large scales (Van Waerbeke et al. 2001, 2002; 
Hoekstra et al. 2002; Bacon et al. 2002; Refregier et al. 2002, and references therein). 

More recently, using new X-ray cluster samples and different virial mass versus tempera- 
ture relations (which are critical for a precise determination of the cosmological constraints), 
cluster normalizations that are either considerably lower (by ~ 2cr; ag^m'^'^ — 0.4 ± 10% 
) or higher (~ 0.6 ± 10%) than the above value have been reported (Borgani et al. 2001; 
Ikebe et al. 2002; Seljak 2002; Reiprich & Bohringer 2002; Viana, Nichol, & Liddle 2002; 
Pierpaoh et al. 2001; here we converted all the relations to the same power-law slope of 0.5 
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for easier comparison). An accurate determination of this parameter is important for two 
reasons. First, the normahzation cTg^ enters exponentially in the evolution of structure in 
the universe; a 20% change in cTg has a significant (exponential) effect on the evolution of 
structure with time and, of course, on the amount of bias in the universe (i.e., how mass 
traces light). Second, if we know as or Qm from other observations, the above relation can 
be used to determine the second parameter. For example, if erg ~ 1, as suggested by some 
observations, then the implied value of differs by nearly a factor of two depending on 
whether the cluster normalization relation is crgfim^'^ — 0.6 ± 10% or 0.4 ± 10%; these values 
imply Q„= 0.36 ± 20% or 0.16 ± 20%, respectively. 

Most of the previous analyses, which use the cluster temperature function, employ a 
smaller number of clusters, and assume a relation between virial cluster mass and temper- 
ature, which sensitively affects the results. The cosmological interpretations are generally 
based on comparisons with theoretical approximations such as the Press- Schechter formal- 
ism. 

In this paper we use the early commissioning data from the Sloan Digital Sky Survey 
(SDSS: York et al. 2000; Stoughton, et al. 2002) to determine a preliminary mass function of 
nearby clusters of galaxies and derive its cosmological constraints. The data cover about 400 
dcg^ with ~300 clusters at z = 0.1 - 0.2 from each of two independent samples (~600 clusters 
in total) — considerably larger than previous samples. The analysis does not use cluster 
virial masses, nor the virial mass temperature relation, which are more difficult to determine 
obscrvationally. Rather, we use cluster masses observed within a fixed radius, as calibrated 
from the observed cluster luminosities and tested against cluster velocity dispersion. We 
use two independently selected cluster samples, identified by different algorithms — the 
Matched Filter method and the color-magnitude maxBCG method; we find consistent results 
for the two samples. We compare the results directly with large scale (Gpc^) cosmological 
simulations as well as with the Press- Schechter formalism to determine the cosmological 
constraints. 

Finally, we note that the current results are based on a very small fraction (4%) of the 
ultimate SDSS 10^ deg^ survey, which will yield thousands of clusters, many with velocity 

dispersions and weak gravitational Icnsing masses. The present mass function is therefore 
preliminary, intended to show the feasibility of using clusters from SDSS by utilizing the 
early commissioning data; larger and more accurate data will become available from SDSS 
in the near future. 



-4- 



2. Cluster Selection from SDSS Commissioning Data 

The SDSS (York et al. 2000) is a 5-band CCD imaging survey that will cover, when 
complete, 10^ deg^ of the high latitude North Galactic Cap, and a smaller deeper region in 
the South, followed by an extensive multi-fiber spectroscopic survey. The imaging survey is 
carried out in drift-scan mode in five SDSS filters, u, g, r, i, z, to a limiting point source 
magnitude of r<23 (Fukugita et al. 1996; Gunn et al. 1998; Lupton, et al. 2001; Hogg et al. 
2001; Smith et al. 2002; Pier et al. 2002). The spectroscopic survey will target nearly one 
million galaxies to approximately r<17.7, with a median redshift of z~0.1 (Strauss, et al. 
2002), and a smaller deeper sample of ~ 10^ Luminous Red Galaxies to r~19 and z~0.5 
(Eisenstein, et al. 2001). 

In this paper we use 379 deg^ of the early commissioning data of SDSS imaging, covering 

the area a(2000) = 355.0°to 56.0°, 5(2000) = -1.25°to 1.25°; and a(2000) = 145.3°to 236.0°, 
5(2000)= -1.25°to 1.25° (runs 94/125 and 752/756; Stoughton, et al. 2002). Clusters of 
galaxies were selected from these imaging data using, among others, a Matched-Filter method 
(Kim et al. 2002a,b) and an independent color-magnitude maximum-likelihood Brightest 
Cluster Galaxy method (maxBCG; Annis et al. 2002). These methods are briefly described 
below. A detailed comparison between these independent cluster selection methods and 
their properties is given in Bahcall et al. (2002). Here we use clusters selected from these 
techniques to determine a preliminary mass function of nearby clusters of galaxies. 

The Matched Filter method HMF (Hybrid Matched Filter; Kim et al. 2002a) is a hybrid 
of the Matched Filter (Postman et al. 1996) and the Adaptive Matched Filter techniques 
(Kepner et al. 1999). This method identifies clusters in imaging data by finding peaks in 
a cluster likelihood map generated by convolving the galaxy survey with a filter based on 
a model of the cluster and field galaxy distribution. The cluster filter is composed of a 
projected density profile model for the galaxy distribution (Plummer law profile), and a 
luminosity function filter (Schechter function), using the typical parameters observed for 
galaxy clusters (within a radius of 1 hr^ Mpc). The HMF method identifies the highest 
likelihood clusters in the imaging data and determines their best-fit estimated redshift (zesj) 
and richness (A); the best-fit richness is proportional to the total cluster maximum hkelihood 
luminosity within a radius of 1 Mpc. A relatively high threshold has been applied to 
the HMF cluster selection [a >5.2, Kim et al. 2002a); therefore, the selected clusters have 
typical richness of A> 20-30 (i.e., h^i{< Ih'^ Mpc) > 20L* ~2xl0^^ h''^ Lq). This threshold 
corresponds to clusters poorer than Abell richness class 0. (For more details see Kim et al. 
2002a). 

The maxBCG method (Annis et al. 2002) is based on the fact that the brightest cluster 
galaxy (BCG) generally lies in a narrowly defined space in luminosity and color (see, e.g. 
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Hoessel & Schneider 1985; Gladders & Yee 2000). For each SDSS galaxy, a BCG hkehhood 
is calculated as a function of redshift based on the galaxy color {g-r and r-i) and magnitude. 
The cluster likelihood is then weighted by the number of nearby red galaxies (located within 
1 Mpc projected separation) that are within the color-magnitude region expected for 
the relevant cluster E/SO galaxy ridgeline. This combined likelihood is used for cluster 
identification. The likelihood is calculated as a function of redshift from z = to 0.5, at 0.01 
intervals. The best estimated redshift is that which maximizes the cluster likehhood. Since 
elliptical galaxies possess very regular colors, they provide excellent photometric redshift 
estimates for their parent clusters. The richness estimator, Ngai, is defined as the number of 
red E/SO ridgehne member galaxies that are brighter than Mj = -20.25 (i.e., 1 mag fainter 
than L*; h = 1), and are located within a 1 Mpc radius of the BCG. (For more details 
see Annis et al. 2002). 

The HMF cluster catalog contains clusters with richness A> 20 and redshift z^st < 0.5 
(Kim ct al. 2002a,b). The selection function for this sample has been determined using 
simulated clusters (see above references). The HMF redshift uncertainty is determined to 
be az = 0.03 (by comparison with measured rcdshifts, Bahcall et al. 2002); the redshift 
uncertainty of the maxBCG clusters is = 0.02. 

In this paper we determine the abundance of HMF clusters as a function of richness for 
nearby clusters (z = 0.1 - 0.2) and use the observed richness - mass relation to determine 
a preliminary mass function for the HMF clusters. A similar analysis is carried out for the 
independently selected maxBCG clusters and the results compared. Each of the independent 
samples contains ~300 clusters within the redshift (z = 0.1 - 0.2) and richness (A> 30 and 
Ngai> 10) hmits used in this analysis. 

3. The Cluster Mass Function 

3.1. HMF Cluster Mass Function 

We determine the mass function of nearby clusters of galaxies using HMF clusters with 
richness A> 30 and redshift z = 0.1 - 0.2. (At z < 0.1, the number of clusters is small and their 
selection less effective; we thus restrict our analysis to the above range.) To minimize false- 
positive detections we use the VCl sample (Visually Confirmed sample, Kim et al. 2002a,b) 
which contains > 80% of all A> 30 HMF clusters, increasing to > 90% for A> 50 clusters. 
The total number of VCl clusters observed within this redshift and richness range is 294 
(uncorrected for selection function). Each cluster is corrected by the appropriate selection 
function for the given cluster richness and redshift as determined from cluster simulations 
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(Kim et al. 2002a). The cluster abundance as a function of richness, from A> 30 to A> 
70, is obtained by dividing the above volume-hmited corrected cluster count by the relevant 
volume (z = 0.1 - 0.2). A fiat Qm= 0.3 cosmology is assumed for the volume calculation, and 
a Hubble constant of Hq = 100 h km/s/Mpc is used. (When fitting to different cosmologies 
in Section 4, the proper self-consistent cosmological volume is used for each flrn value.) 

Two corrections arc applied to the cumulative cluster richness function. First, wc correct 
the abundance of clusters above a given richness, n (> A), for the effect of redshift uncertainty 
in the HMF clusters, = 0.03 (see Section 2). The correction factor is determined using 
Monte Carlo simulations of realistic cluster distribution with redshift and richness, which 
is convolved with the observed Gaussian scatter in redshift, az = 0.03. We find that the 
redshift uncertainty has a small effect, causing about 10% more clusters to be scattered into 
the z = 0.1 - 0.2 volume than are scattered out. We correct the cluster abundances downward 
by this small correction. Second, we correct the derived cluster abundance for the effect of 
uncertainty in the HMF richness, estimated to be 20% based on cluster simulations. We 
use Monte Carlo simulations with a realistic richness function, convolve it with the known 
observational selection function to produce the observed number of clusters as a function of 
true richness and then scatter the richness with the observed uncertainty to yield the observed 
richness function. Comparing the observed and true richness functions in 10^ simulations 
we determine the proper correction factors and their dispersion, which we apply to the data. 
We find that the observed abundances are larger than the true ones, as expected due to the 
excess scatter of the more numerous low richness clusters to higher richness; this effect is 10% 
at A~ 30 - 40, increasing to 35% - 55% at A~ 60 - 70. We correct the cluster abundances 
for this effect, and use the observed variance in the final error analysis discussed below. 

The uncertainties in the observed cluster abundance include the statistical uncertainties 
(Na), the uncertainties in the selection function (±15%) and in the false-positive correction 
(±15%), and the uncertainties derived from the Monte Carlo simulations for each of the two 
corrections above (the redshift correction factor has an uncertainty of 4% to 42% for the 
range A>30 to A>70, and the richness correction factor has an uncertainty of 3% to 23% 
for the same range). 

To determine a cluster mass function from the above cluster richness function we need 
to convert the cluster richness thresholds to a mean cluster mass. Throughout this paper we 
use cluster mass within a given fixed radius (not virial mass); this mass is more accurately 
obtained from observations since the virial radius is not precisely known. We convert richness 
to mass in two independent ways, both from observations. First, we use the mean cluster 
luminosity measured in the SDSS data for all clusters stacked as a function of their richness. 
The cluster luminosity is observed within a radius of 0.6 h"^ Mpc, in the r-band, for galaxies 
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brighter than = -19.8 (K-corrected for each galaxy type following Fukugita et al. 1996), 
and corrected for a similarly determined local background in five separate locations (which 
allows us to account for the variance in the background correction; Hansen et al. 2002; 
Bahcall et al. 2002). We use the mean observed cluster luminosity Lo.e for clusters with 
richness threshold of A=30, 40, 45, 50, 60 and 70. The observed mean luminosity Lo.e of 
the stacked clusters is presented as a function of richness in Figure 1. We note that any 
biases or uncertainties in the richness parameter (e.g., Kim et al. 2002a) are calibrated 
out in this procedure since the actual mean cluster luminosities are directly measured by 
this method. The richness parameter serves only as a tracer; a richness bias will properly 
calibrate itself by the measured mean luminosity (as is in fact seen by the non-linear relation 
between Lg.e and A). The cluster luminosity is corrected to include the unobserved faint-end 
of the cluster luminosity function, that is, all galaxies fainter than M^. = -19.8. For the 
observed SDSS Schechter luminosity function parameters of the HMF clusters (within 0.6 
h-^ Mpc), a = -1.08 ± 0.01 and M; = - 21.1 ± 0.02 (Hansen et al. 2002; see also Goto 
et al. 2002; h = 1), we adopt a correction factor of 1.42 ± 0.08 for the added contribution 
of faint galaxies to the total HMF cluster luminosity. The cluster mean luminosity is then 
converted to cluster mass, M(< 0.6 Mpc physical), using the mean observed cluster 
M/Lr ratio for each richness threshold (Bahcall & Comerford 2002). The observed best- 
fit M/L is used (based on the means of 20 clusters and 33 groups): M/L^^tot (z=0) = 
142±32 + (23±5) Tkev h (Bahcall & Comerford 2002). The mild increase of M/L with 
temperature, T^g^,, (seen both in observations and in simulations, e.g., Bahcall et al. 2000), is 
accounted for at each richness threshold using the observed correlation between richness and 
velocity dispersion (see below) and the observed mean relation between velocity dispersion 
and temparature (cr^= 332 T^-^^kms^^, Lubin & Bahcall 1993). This effect is small for the 
range of cluster temperatures studied here (T ~ 0.9 to 4 kev). The mean observed M/L^^tot 
at z—0 is converted to M/ Lj. (where is the relevant SDSS Petrosian r luminosity) using 
the conversions given by Fukugita et al. (1996), Bahcall & Comerford (2002), Strauss, et 
al. (2002). We use Lr = 0.85 Lr,tot (Strauss, et al. 2002, for ~ 60 - 70% of cluster light 
contributed by early type galaxies), M/Lr^tot = 0.94 M/L^^tot, and M/Lj. (z = 0.17) = 0.943 
M/Lr (z = 0) (Carlberg et al. 1997a; Bahcall & Comerford 2002). The above yields M/Lr 
values (at z=0.17) that range from 170 at A>30 to 235 at A>70. The mean mass of clusters 
(within 0.6 h^^ Mpc) is then determined for the relevant richness thresholds (A>30 to A>70). 

The uncertainties in the mean mass estimates are derived from the combined uncer- 
tainties in the observed mean luminosity-richness relation, the uncertainty in the observed 
mean M/ L ratio, and the smaller uncertainties in the corrections applied above. In order to 
determine the proper uncertainties in the luminosity-richness relation that are relevant for 
the mean cluster mass estimates, we generate 10^ Monte Carlo simulations with a realistic 
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L-A relation and richness function. We introduce a Gaussian redshift scatter of cr^ = 0.03 as 
well as a ~ 20% to 30% uncertainty in individual cluster luminosities. We recover the mean 
"observed" L-A relation from the lO"' simulations and the mean luminosities at the relevant 
cluster richness thresholds (for the observed redshift range). The recovered mean relation is 
consistent with the input L-A relation. The derived Icr variance in the recovered L-A relation 
from the simulations ranges from 11% at A=30 to 25% at A=70. We use these uncertainties 
in estimating cluster mass uncertainties. The uncertainty in the mean observed M/ L ratio, 
15% for the relevant cluster richnesses, is combined with an additional 8% uncertainty in 
the conversion factors described above and a 6% uncertainty in the luminosity function faint 
end extrapolation (see above). The mass uncertainties thus range from 20% at A~30 to 31% 
at A~70. The cluster abundances have been corrected for this scatter using Monte Carlo 
simulations, as described above for the richness-function abundance correction. 

For comparison with other commonly used cluster masses, as well as for direct com- 
parison with available cosmological simulations, we also determine the mass function for 
cluster masses within two additional frequently used radii: the slightly smaller radius of 0.5 

Mpc, using the observed mean luminosities L0.5, and, for illustration purposes, also the 
larger comoving radius of 1.5 h^^ Mpc; the latter is obtained by extrapolating the 0.6 
Mpc luminosity to 1.5 h^^ Mpc comoving radius (= 1.28 Mpc at z = 0.17) using the 
typical observed luminosity profile in clusters { Pl ^ for R < R200 and ~ R~^-^ for R > 
R200, where R200 is the radius within which the cluster overdensity is 200 times the critical 
density; Carlberg et al. 1997b; Fischer & Tyson 1997). 

These two mass functions are compared in Figure 2 with the mass function obtained 
from large scale cosmological simulation (Bode et al. 2001) of the concordance LCDM model 
(Bahcall et al. 1999): 0.3, A= 0.7, crg^ 0.9 (i.e. crgQ^n"-^ = 0.49), and /i= 0.67. This 

simulation used a 1 Gpc box size and 1024^ dark matter particles, with a particle mass 
of 2.3 X 10^° h-^ Mq, and softening length of 14 h'^ Kpc (Bode et al. 2001). Such a large 
box ensures a statistically valid sample of simulated clusters, and the high particle number 
ensures that the clusters are well resolved — for the smallest clusters considered here there 
are over 10'^ particles within 0.5 h^^ Mpc. The details of the simulation and the method 
of computing the mass function are described in Bode ct al. (2001). The simulated mass 
function is presented as a function of M(< 0.5 h'^ Mpc physical) and M(< 1.5 Mpc 
comoving) at z = 0.17, for direct comparison with the observations. Figure 2 shows that 
the shape of the SDSS mass function agrees well with that expected from the cosmological 
simulations but the normalization is significantly lower than expected from the concordance 
model. The best-fit function, with a lower as-^m amplitude, is also presented in Figure 2; 
it is discussed in Section 4. 
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The observed HMF cluster mass function for M(<0.6 Mpc) is presented in Figure 
3. As a further consistency test, we estimate mean cluster masses using an entirely indepen- 
dent method: the observed correlation between mean cluster richness and cluster velocity 
dispersion. We use cluster velocity dispersions of 19 clusters determined from the SDSS 
spectroscopic survey (for clusters with ~30 to 160 redshifts) as well as from several Abell 
clusters available in the literature (Mazure et al. 1996; Slinglend et al. 1998; Abell 168, 295, 
957, 1238, 1367, 2644). Even though the number of clusters with measured velocity disper- 
sion is not large and the scatter considerable, a clear correlation between median velocity 
dispersion and richness is observed, as expected; we find a best-fit relation av{km s~^) ~ 
10.2 A. We estimate mean cluster mass (within 0.6 Mpc) from this relation and use it 
to illustrate consistency with the mass function determined from the entirely independent 
cluster luminosity method discussed above. We use the observed relation between cluster 
mass and cluster velocity dispersion derived from observations of weak gravitational lensing 
of clusters: M(< 0.6 Mpc) = 0.0717 ks aji^f 10^^ h'^ M© (where aioo is in 100 km s'^ 
Hjorth et al. 1998; also Bahcall & Sctte 2002). This relation is obtained from the observed 
relations M/R — 0.88 ks T(kev) for R < 1 Mpc, where ks is the small overdensity cor- 
rection factor {ks — 0.76, 0.9, 1, 1.1, 1.15, respectively, for cluster overdensity oi 6 — 100, 
250, 500, 1000, 2500; see references above and Evrard et al. 1996), and a-t,{kms~^) — 332 
T^g^ (Lubin & Bahcall 1993). The cluster mass function determined from this independent 
method, performed as a consistency check, is in full agreement with the mass function deter- 
mined earlier using cluster luminosities; the results are compared in Figure 3. The velocity 
dispersion comparisons from the two methods — i.e., the velocities inferred from the clus- 
ter luminosity-mass method and the directly observed velocity dispersions are shown as a 
function of richness in Figure 4. The excellent agreement between these two independent 
methods supports the mass determination discussed above. 



3.2. maxBCG Cluster Mciss Function 

For comparison, we also determine the cluster mass function from the independently 
selected maxBCG clusters. This method uses a completely independent selection criterion: 
the maxBCG selection technique assumes no cluster filters or profiles; rather, it selects 
clusters based on the red colors and magnitudes of the brightest cluster galaxies (Section 2). 
A comparison of the two mass functions can therefore provide further support for the above 
results. 

We follow the same procedure for the maxBCG clusters as described above for the 
HMF clusters. We use the observed mean luminosity Lo.e of all stacked maxBCG clusters 
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as a function of richness, Ngai (where Ngai is the maxBCG cluster richness, Section 2); 
the data are presented in Figure 1. Wc extrapolate the luminosity to the faint-end of the 
cluster luminosity function (within 0.6 h'^ Mpc; a = -1.05 ± O.OL M* = -21.25 ± 0.02 for 
the maxBCG clusters for h = 1, Hansen et al. 2002; see also Goto ct al. 2002), yielding a 
correction factor of 1.34 ± 0.06, and convert the cluster luminosity to mean cluster mass 
using the mean observed M/L ratios. All maxBCG clusters (357 clusters) with richness 
Ngai> 10 (comparable in richness to HMF clusters with richness A> 30; Bahcall et al. 2002) 
in the redshift range z = 0.1 - 0.2 are used. Corrections and uncertainties are calculated as 
described above (with az— 0.02, ANgai— 10% - 15%). The selection function and the false- 
positive correction factor for the Ngai>10 maxBCG clusters at z = 0.1 - 0.2 are estimated 
from simulations to be ~ 0.9 - 1 ± 15% each. The fraction of HMF clusters that are 
found by the maxBCG method is 61% (for maxBCG matches with Ngai>Q located within 
1 h^^ Mpc projected separation). This is consistent with the maxBCG selection function 
(~85% for Ngai'>6), the HMF false-positive detection rate (~20%, for A>30), and the smaller 
correction due to redshift uncertainty (see Bahcall et al. 2002). The overlap rate decreases 
considerably if only Ngai>10 maxBCG matches are considered for the HMF clusters (within 
1 separation and ±0.05 in redshift). This is as expected due to the large scatter of 
clusters across the richness threshold; Monte Carlo simulations indicate that a richness cut 
reduces the overlap by ~ 55% - 60% (Bahcall et al. 2002). 

There are 58 Abell clusters located within the current survey area; all 58 clusters are 
detected by the combined HMF and maxBCG clusters (maxBCG finds all 58 clusters; HMF 
finds 49 of the 58 clusters, consistent with the respective selection functions). A few of the 
clusters are detected below the thresholds used here (i.e., A<30, Ngai<lO, z<0.1, or z>0.2). 
In addition, there arc 5 X-ray clusters from the XBACs sample (Ebcling ct al. 1996) located 
within the survey region. All 5 are detected by the HMF and maxBCG methods. 

The maxBCG cluster mass function is presented in Figure 5; it is superposed, for com- 
parison, on the HMF mass function from Figure 3. A good agreement between the two 
independent mass functions is observed. This agreement provides further support to the 
above estimate of the SDSS cluster mass function. In addition, we use the best-fit rela- 
tion observed between mean cluster velocity dispersion and richness for maxBCG clusters 
(21 clusters with measured dispersions) as an additional test; the observed median relation, 
(T^(km/s) = 93 Ngai^'^^ (Figure 4), is used to independently estimate cluster masses using M 
(< 0.6 h-^ Mpc) = 0.0717 ks al^J 10^^ h'^ M©, as described above. The results of the two 
methods (Figure 5) are consistent with each other. 

We can compare the observed cluster mass function with the recently observed cluster 
temperature function by Ikebe et al. (2002). For this purpose we use the observed relation 
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discussed above, M (< 0.6 h"^ Mpc) = 0.53 ks T(kev) 10-^^ M©, to derive an approximate 
temperature function from the above mass function. We find a good agreement between the 
mass function (based on 300 clusters at z = 0.1 - 0.2) and the Ikebc ct ah temperature 
function (based on 60 X-ray clusters, mostly at z < 0.1) (accounting for the slightly higher 
redshift of the SDSS sample, in accord with cosmological simulations; Bode et al. 2001). The 
SDSS mass function reaches to poorer clusters, of lower temperature (T ~ 1 kev), as com- 
pared with the X-ray temperature function (T > 2 kev); because of the small area covered, 
the current SDSS sample does not contain the most massive clusters — these highest mass 
clusters will become available as the sample size increases. The agreement between these 
independent determinations provides further support of the current cluster mass function 
results. 



4. Cosmological Implications 

The cluster mass function places one of the most powerful constraints on the cosmolog- 
ical parameters and as; it determines the important cluster normalization relation, i.e., 
the value of as as a function of Q^. 

Early data of the cluster mass function (Bahcall & Cen 1992), and the cluster tem- 
perature function (Henry & Arnaud 1991; White, Efstathiou, & Prenk 1993; Eke, Cole, & 
Prenk 1996; Eke et al. 1998; Pen 1998), provided a cluster normalization relation of crsilm^'^ 
~ 0.5 (± ~10%; see above references for details). This powerful relation implies that for 
fT8~ 1, 0.25 (with slight differences depending on a flat versus open cosmology). For 

Qrn= 1) the required normalization of as= 0.5 implies a strong bias which is not supported 
by observations (Bahcall, Lubin, & Dorman 1995; Bahcall et al. 2000; Bahcall & Fan 1998; 
Feldman et al. 2001; Verde ct al. 2002; Lahav et al. 2002). 

More recently, using different X-ray cluster samples and different relations between virial 
mass and cluster temperature, a somewhat lower normalization value has been suggested 
(Borgani ct al. 2001; Ikcbc ct al. 2002; Seljak 2002; Reiprich & Bohringer 2002; Viana, 
Nichol, & Liddlc 2002), although higher values have also been reported (Pierpaoli et al. 
2001). In this section we compare the preliminary mass function of SDSS clusters with 
analytic predictions to determine the best-fit cosmological parameters. 

The mass function for a given cosmology can be predicted using the analytic formalism 
of Press & Schechter (1974), as in for example Eke, Cole, & Frenk (1996); Kitayama & Suto 
(1996); Viana & Liddle (1996); Henry (2000); and Reiprich & Bohringer (2002). While 
fairly successful in matching the results of N-body simulations, the standard P-S formahsm 
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tends to predict too many low mass clusters and too few higher mass clusters. An improved 
fitting formula which better reproduces the results of N-body simulations is given by Jenkins 
et al. (2001); we will use this in preference to the standard P-S formula. 

For a given set of cosmological parameters, we begin by calculating the linear matter 
power spectrum using the publicly available CMBFAST code (Seljak & Zaldarriaga 1996). 

Knowing the power spectrum, we calculate the variance of the linear density field and thus 
find the mass function with equation B3 of Jenkins ct al. (2001). As Hu & Kravtsov (2002) 
show, this formula is appropriate for a definition of cluster mass within a sphere enclosing 
a mean ovcrdcnsity which is a fixed multiple of the mean density. Because the HMF and 
maxBCG masses are instead derived within a fixed radius (independent of density or mass), 
we adjust the analytic masses using the mass distribution corresponding to the NFW density 
profile (Navarro, Frenk, & White 1997), which provides an accurate representation of N-body 
results. To do this we follow the method presented in the Appendix of Hu & Kravtsov (2002). 

The resulting analytic prediction, using cluster mass within two different radii of 0.5 
Mpc (physical) and 1.5 Mpc (comoving), is shown for the concordance model in Figure 
2 as a dotted line. This analytic function is in excellent agreement with the direct N-body 
simulation (dashed line) over the relevant range of masses (though at higher or lower masses 
the agreement may not be so close). 

To determine the best-fit mass function and the implied values of and as, we compare 
the differential binned mass function to the above theoretical prediction using a standard 
procedure. The last data point in each sample is not included in the best fit determinations 
since this point contains only a few clusters (~7) with considerably less well determined mass. 
Only spatially fiat models are considered. The Hubble constant is kept fixed at Hq = 72 km 
s^^ Mpc^-*^ (Freedman et al. 2001), with a baryon density of Q^h'^ = 0.02 (Buries, Nollctt, 
& Turner 2001), and CMB temperature Tcmb = 2.726 (Mather et al. 1994). We assume a 
primordial power spectral index n = 1. 

The results, presented in Figures 2, 3 and 5, show a good agreement between the 
shape of the observed and theoretical LCDM mass function, but the observed function 
has a significantly lower normalization than the canonical value of as^rn''^ = 0.5 (the latter 
indicated by the dashed and dotted curves in Figure 1, representing the concordance LCDM 
model with ftm— 0.3 and as— 0.9). Model simulations with a somewhat lower value of 
asO,jn^-^ — 0.45 (for a flat Quintessence model and Open CDM, see Bode et al. 2001; not 
shown here to avoid crowding) also exhibit higher normalization than observed. 

The best-flt mass function is presented in Figures 2 and 3 for the HMF clusters, and in 
Figure 5 for both the maxBCG and HMF clusters; the two independent best flt functions yield 
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similar results (Figure 5). The cosmological constraints derived from the minimization are 
summarized in Figure 6, showing the allowed parameter range of l^m-o"8 for both the HMF 
and the maxBCG samples. The best-fit parameters are ^ra= 0.175 ±007, 0^= 0.92 ±020 
(l-d) for the HMF clusters, and 0.195 ±[!:[}?, ^8= 0.9 ±[]J for the maxBCG clusters 

(Figures 2, 3 and 5). The best-fit contours in Figure 6 show that high values of fim.(> 0.4) 
and low values of (y%{< 0.6) are ruled out by the data at > 2-(t; these yield mass functions 
that are too steep compared to the data. On the other hand, low values of ^^(down to ~ 
0.1) and high values of erg (up to ~ 1.2) are supported by the data. These results are obtained 
using the M(< 0.6 hT^ Mpc) mass functions; similar results are obtained with M(< 1.5/i~^ 
Mpc comoving). 

The mean best-fit parameters of the observed mass function (with \-a) are: 

= 0.19±n? (1) 

(78 = 0.9±°i (2) 

The best-fit normalization relation (Figure 6) is: 

(78^^^ = 0.33 ±0.03 (0.1<Q^<0.4) (3) 



or, equivalent ly. 



= i^-^r ± 10% . (4) 



For comparison with previous results, this relation (3) is 20% lower than the standard 
normalization value of as^^^-^ — 0.5 (or similarly, as^m^'^ — 0.44). As discussed earlier, 
this conclusion has a non-negligible implication for as and ilrn as seen in equations (1) and 
(2) for the best-fit parameters (and discussed below). 

The results are consistent with the recent temperature function results of Ikebe et al. 
(2002) and Seljak (2002). The higher amplitude obtained by some of the earher work is most 
likely due to a combination of factors including use of the uncertain and sensitive theoretical 
(rather than observational) mass-temperature relation (when applied to X-ray clusters), 
smaller sample size, and overestimated cluster abundance (in some optical samples). Recent 
weak Icnsing observations on large scales yield results that range from (Ts^m"'^ ~ 0.4 to 0.6 
(or agflm'^'^ ~ 0.34 to 0.52 when converted to a power law slope of 0.6, at flm^ 0.25, for 
easier comparison; see references in Secrtion 1); the low end of this range is consistent with 
the current results, but the high normalization values reported are inconsistent at the 2-a 
level. A complementary analysis of the maxBCG cluster halo occupation function (Annis et 
al., in preparation) yields consistent results with those obtained above. 
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Our best-fit as-^m constraints are compared with previous results in Figure 7. This 
comparison illustrates the agreement of recent temperature function results with the current 
constraints, and shows the wide range among the earlier, higher as{^m) normalization results. 
The weak lensing analyses yield constraints that lie mostly at the higher (Ts{Qrn) range; the 
lensing results of Hoekstra et al. (2002) are consistent with our current constraints. 

This new cluster normalization has important implications for Q^. It is frequently 
assumed that is 0.3, and the above relation is used to determine cig (thus typically 
referred to as "low normalization", implying a lower than expected (jg value for ilm= 0.3). 
However, the value of Qm is not accurately known, and could be as low as ~ 0.15 (see, e.g., 
Bahcall et al. 2000). At the same time, there are several measurements that suggest that 
(78~ 0.9 - 1. This "high" normalization is obtained from the very mild evolution observed 
in the cluster abundance to z~l (Bahcall & Fan 1998; Donahue & Voit 2000); the flattening 
of the M/L function on large scales (Bahcall, Lubin, & Dorman 1995; Bahcall et al. 2000); 
SDSS observations of the galaxy power spectrum on large scales (Szalay et al. 2002); and 
the observation of no-bias in the galaxy distribution in the 2dF and infrared surveys (Verde 
et al. 2002; Feldman et al. 2001; but see, however, Lahav et al. 2002). This observationally 
suggested normalization of (Tg— 0.9 - 1 is fully consistent with our best-fit value for the SDSS 
cluster mass function (equations 1 and 2); it implies a low mass density of flm— 0.19. This 
is also consistent with the low value indicated by the M/ L function on large scales. 



5. Conclusions 

We determine the mass function of nearby clusters of galaxies using ~ 300 clusters at z 
= 0.1 - 0.2 selected from 400 deg^ of early SDSS commissioning data. Two independent 
cluster samples are used based on the Matched Filter and the color-magnitude maxBCG 
methods. The two samples yield consistent results. The analysis uses cluster masses within 
a fixed radius. The mass function is compared directly with large, Gpc^ cosmological sim- 
ulations. We find a best-fit cluster normalization relation of (TgQ^^-^ = 0.33 ± 0.03, or 
equivalently ag— 10%. This result is significantly lower than the previous canoni- 

cal value of CTgQ^^'^ = 0.5 ± 0.05. The shape of the cluster mass function partially breaks 
this degeneracy; we find best-fit parameters of Qm= 0.19 ±007 and 0"g= 0.9 ±02- These 
values are consistent with the independent observationally suggested normalization of (Tg~ 
0.9 - 1 observed from cluster abundance evolution (Bahcall & Fan 1998; Donahue & Voit 
2000), the flattening of the M/L function on large scales (Bahcall et al. 2000), and SDSS, 
2dF, and infrared large-scale structure observations (Szalay et al. 2002; Verde et al. 2002; 
Feldman et al. 2001). 
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These preliminary results from early commissioning data of 4% of the ultimate SDSS 
survey will be greatly improved as additional SDSS data become available for thousands of 
clusters, many with measured velocity dispersions and weak lensing masses. Improvements 
in the sample size and accuracy, and in the relevant scaling relations are needed in order 
to achieve greater precision in the determination of the mass function and its cosmological 
implications. 

The Sloan Digital Sky Survey is a joint project of The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins Uni- 
versity, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), 
the Max-Planck- Institute for Astrophysics (MPA), New Mexico State University, Princeton 
University, the United States Naval Observatory, and the University of Washington. Apache 
Point Observatory, site of the SDSS telescopes, is operated by the Astrophysical Research 
Consortium (ARC). Funding for the project has been provided by the Alfred P. Sloan Foun- 
dation, the SDSS member institutions, the National Aeronautics and Space Administration, 
the National Science Foundation, the Department of Energy, the Japanese Monbukagakusho, 
and the Max Planck Society. The SDSS Web site is http:/ /www. sdss.org. 
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Fig. 1. — Observed cluster luminosity versus richness for HMF and maxBCG clusters. Clus- 
ter luminosity is observed in the r-band, within a radius of 0.6 Mpc, for stacked clusters 
at a given richness. The luminosities are k-corrected, background subtracted, and integrated 
down to Mr — -19.8. Dark squares represent binned data. (See section 3) 
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Fig. 2. — The HMF cluster mass function for masses within radii of 0.5 Mpc (left panel) 
and 1.5 h^^ Mpc comoving (right panel). The solid line is the best-fit analytic mass function 
(determined at 0.6/i^^ Mpc and extrapolated to the appropriate radius assuming an NEW 
profile), with fi^=0.175 and (78=0.92. In each panel the dashed fine is the mass function 
measured from an N-body simulation of the concordance LCDM model with Jl^=0.30 and 
(78=0.90; the dotted line is the analytic prediction for this simulated cosmology. 
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Fig. 3. — The HMF cluster mass function, showing masses (within 0.6 Mpc) determined 
from both the luminosity - mass calibration (filled circles) and the independent velocity 
dispersion - mass relation (open circles). (The observed cluster abundances assume a volume 
corresponding to a fiat Qm= 0.2 cosmology.) The best-fit analytic model, with flm= 0.175 
and (78= 0.92, is shown by the solid line. 
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Fig. 4. — Relation between observed cluster velocity dispersion a and cluster richness (tri- 
angles are SDSS observed velocity dispersions, circles are Abell clusters, dark squares are 
medians, and solid line is best fit to the velocity data. Stars represent SDSS observation 
of stacked cluster data, shown for comparison only). The median observed velocities are 
compared with the velocity determined from the cluster masses derived from the mean 
richness-luminosity-mass relation used in Section 3 (represented by open squares and the 
best-fit dashed line). A good agreement between the two methods is seen. 
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Fig. 5. — The maxBCG and the HMF cluster mass functions, showing masses determined 
from both luminosity - mass relation (solid triangles: maxBCG; solid circles: HMF) and 
velocity dispersion - mass relation (open triangles: maxBCG; open circles: HMF). The best- 
fit analytic models are shown by the dashed line (maxBCG; fljn= 0.195, as= 0.90), and solid 
line (HMF; f]„= 0.175, a8= 0.92; as in Figure 2). 
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Fig. 6. — Allowed Vtrn^a^ range: one- and two-cr confidence contours for HMF (solid lines) 
and maxBCG clusters (dashed lines). The dotted curve is the best-fit relation erg = 0.33fi~°-^ 
— (fr^)°'^- '^^^ best-fit Q,rm CTs values are shown by the dark circle (HMF) and triangle 
(maxBCG). 
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Fig. 7. — Comparison with previous results: the allowed ^m-'^s range from our present SDSS 
mass function is compared with previous results from cluster temperature functions. The 
current 2-a constraints (from Fig. 6) are represented by the shaded area, with the best-fit 
relation shown by the solid fine, and the best-fit parameters {0,^— 0.19, as— 0.9) indicated 
by the filled circle. The various dotted and dashed curves represent the best-fit relations 
given by previous work (for clarity, the allowed width of each range is not shown; it typically 
corresponds to ± ~10% in cig at a given Qm)- Dotted lines represent early results (White, 
Efstathiou, & Frenk 1993; Viana & Liddle 1996; Eke, Cole, & Frenk 1996; top to bottom); 
short-dashed lines are mostly re-analyses of same/similar early data (Pierpaoli et al. 2001; 
Oukbir & Arnaud 2001; Proty Wu 2001; top to bottom); long-dashed lines represent lower 
normalization relations obtained from recent analyses and/or recent samples (Seljak 2002; 
Reiprich & Bohringer 2002; Viana, Nichol, & Liddle 2002; top to bottom at higher Q^). 
Weak lensing analyses on large scales (Sections 1, 4; not shown) yield results that are mostly 
near the high range of these curves, with Hoekstra et al. 2002 fitting well on our best-fit 
relation. 



